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With the release of the data from the Boomerang and MAXIMA-1 balloon nights, estimates 
of cosmological parameters based on the Cosmic Microwave Background (CMB) have reached un- 
precedented precision. In this paper I show that it is possible for these estimates to be substantially 
biased by features in the primordial density power spectrum. I construct primordial power spectra 
which mimic to within cosmic variance errors the effect of changing parameters such as the baryon 
density and neutrino mass, meaning that even an ideal measurement would be unable to resolve 
■ the degeneracy. Complementary measurements are necessary to resolve this ambiguity in parameter 

f*"^ ' estimation efforts based on CMB temperature fluctuations alone. 
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I. INTRODUCTION 



With the release of the Boomerang (| and MAXIMA-1 § data sets, the promise of using the Cosmic Microwave 
Background to provide precision constraints on cosmological parameters has become a reality. The inflationary 
prediction of a flat universe has received spectacular confirmation, and remarkably good estimates of other cosmological 
parameters have emerged 1^,0. More detailed observations of the CMB, in particular NASA's MAP satellite || and 
the ESA's Planck Surveyor ||, should allow for exquisitely precise determination of a handful of "fundamental" 
cosmological parameters such as the baryon density and the Hubble constant, 
l/^ ■ However, current parameter estimation efforts based on observation of the CMB suffer from an inherent ambiguity: 
\ lack of knowledge about the primordial density fluctuation spectrum. Primordial density and gravity wave fluctuations 
are the underlying source of the fluctuations in the CMB. The ability to predict the form of the CMB temperature 
anisotropy depends on knowledge of the form of these fluctuations. The standard assumption is that the density 
fluctuation power spectrum is a featureless power law, P (k) oc k n . Such an assumption is theoretically supported 
CLi by inflation, which predicts such a power law spectrum. It is also empirically supported by observation of large 
scale structure. Current large scale structure data, however, only loosely constrain the form of the primordial power 
spectrum on scales relevant for CMB anisotropy measurements, k < 0.1 hMpC^ 1 0. A fully empirical approach to 
the data forces one to consider an arbitrary density power spectrum as input to the parameter estimation problem. 
This in effect adds so many free parameters (and associated degeneracies) that precision determination of parameters 
| from the CMB temperature anisotropy alone becomes an impossibility. In this paper I expicitly construct primordial 
power spectra that mimic the effects of changes in cosmological parameters to within cosmic variance errors, meaning 
that even a perfect measurement of the CMB temperature fluctuations would be unable to resolve the degeneracy. 
This requires large deviations from the standard scale-invariant power law spectrum, 5P(k)/P(k) ~ 0.5, which would 
be detectable with complementary measurements, such as CMB polarization and improved large-scale structure data 
like that from the Sloan Digital Sky Survey M. 



II. THE PRIMORDIAL POWER SPECTRUM AND THE CMB 



All observations of the CMB anisotropy published to date have been maps of the anisotropy in the temperature of 
the CMB, AT /Tq ~ 10~ D , where To = 2.728 K. It is convenient to quantify these temperature fluctuations as a set of 
multipole moments, 

^^=E E a lm Y tm {6A). (1) 
e=om=-e 
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Lack of a preferred direction implies that the amplitudes ae m will be independent of m, 



( a e'm' a em) = CiSu'Smm' , (2) 

where the angle brackets denote an average over realizations. For Gaussian fluctuations, the set of Cts completely 
characterizes the temperature anisotropy. If the fluctuations are non-Gaussian, higher order correlation functions 
are necessary to fully characterize the anisotropy. The power of the CMB in constraining cosmological parameters 
comes from the fact that one is using a large number of parameters (the CV spectrum) to constrain a dozen or so 
"fundamental" parameters such as the baryon density, the ratio of gravity wave perturbations to density perturbations, 
and so on. Any given CMB observation will have access to a range of C/s with error bars determined by a combination 
of the observation's sky coverage, sensitivity and angular resolution. However, there exists a fundamental limit to the 
accuracy with which the Cg spectrum can be measured, referred to as cosmic variance, 



Cosmic variance is simply a finite sample size effect coming from the fact that there is only a single sky to measure, 
and is more significant at long wavelength (smaller £) . Cosmic variance will play a central role in the discussion below 
for the reason that any two Cg spectra which have a small enough \ 2 over cosmic variance errors are in principle 
indistinguishable, even when subject to an ideal measurement. 

Going from a Cg spectrum to cosmological parameters requires assumptions about the form of the density power 
spectrum. A typical assumption when discussing the ability to constrain cosmological parameters using the CMB is 
that the primordial density power spectrum is a power law: 

P(k)ozk n , (4) 

where the spectral index n is close to the scale- invariant value n = 1. The recent Boomerang measurement || yields 
a value n = 0.87+q os f° r a se ^ °f P r i° rs consistent with inflation, and the MAXIMA-1 measurement has a best fit of 
n = 0.99 ±0.09 j|. Such a featureless power law spectrum is in fact predicted by most inflationary models. A natural 
extension of this approximation is to allow "running" of the spectral index, 

'"' *0. (5) 
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Such running is a feature of some inflation models |9|-[l4|| . Future CMB observations will be able to detect running of 
the spectral index as small as dn/dlog(k) ~ 0.05 |f 1| . Taking this a step further, Lesgourges, Prunet and Polarski 



considered CMB constraints on models with broken scale invariance |15|. In principle, however, the primordial power 
spectrum can contain many more free parameters than just one or two. The purpose of this paper is to investigate 
the possibility that features in the primordial power spectrum could mimic the effect of other parameters and thus 
"fool" parameter estimation efforts which assume a featureless power law spectrum. It is to be expected that such 
confusion will be possible, since the number of parameters available in an arbitrary primordial power spectrum is 
much larger than the parameter space considered by typical parameter estimation efforts, and can in principle be 
larger than the number of multipolcs available in the CMB for measurement. I show below that the number of free 
parameters necessary to confuse parameter estimation efforts is fewer than a hundred. In the next section I discuss 
the numerical methods used to approach the problem. 



III. NUMERICAL METHODS 



Boltzmann codes for calculating the CMB spectrum normally take a set of cosmological parameters as input and 
generate a spectrum of Cgs. In this paper I construct primordial power spectra which would cause significant mis- 
estimation of parameters if a power law fluctuation spectrum were to be assumed in the likelihood analysis. This 
amounts to solving the inverse problem - given a particular Cg spectrum, what primordial power spectrum is necessary 
to produce it? I adopt a simple x 2 minimization procedure for solving this problem. The procedure is as follows. 

I choose two models and generate CMB multipole spectra for each using the CMB FAST code |f6| . One I call the 
"true" model, with a set of parameters assumed to represent the actual underlying cosmos. The second I call the 
"target" model, generated with a different set of parameters but the same density power spectrum, which I take to 
be the scale invariant spectrum P (k) cx k. The goal is to deform the density power spectrum in such a way as to 
mimic the Cg spectrum of the target model while retaining the underlying parameters of the true model. The observer 
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will see the target spectrum, but the true parameters of the universe are those of the true model. The input power 
spectrum is binned into Nk = 75 bins i with wavenumber ki. This binning provides a good balance between accuracy 
and computational efficiency. CMBFAST samples the power spectrum at a much higher resolution, so the binned 
data is smoothed using cubic spline interpolation (l7). The "trial" density power spectrum P tlied (k) is initially a 
power law and is adjusted iteratively until the trial CMB spectrum Cf lal matches the target CMB spectrum C e ars ° 
to within cosmic variance errors. The first step is to numerically calculate a correlation matrix between the CYs and 
the primordial power spectrum 



Dm = 



dC) 



'trial 



dP{h) 



(0) 



where ki is the wavenumber of the z'th bin and P (ki) is the amplitude of the density power spectrum at that scale. 1 
Note that this is not a square matrix: I = 1, . . . N( and i — 1, . . . , Nk- I use Nk — 75 and Ng = 1500. It is also 
necessary to compute the derivative of the x 2 between the trial and target spectra: 
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I take the errors to be the cosmic variance limit 
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This is particularly important: any two spectra that agree within cosmic variance limits are impossible to distin- 
guish with any CMB measurement, no matter how sensitive. The result is therefore independent of any particular 
experiment. A new density power spectrum is calculated by gradient descent |l8| : 



P" 



(ki) - F trial (ki) 
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dP(k J 



(9) 



where A is a constant whose value is picked to ensure smooth convergence (I use A = 0.1), and the matrix a^- 
defined in terms of the derivative matrix (^) as: 



is 
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(10) 



Inversion of the matrix otij is in practice somewhat problematic, as the matrix is in numerically singular. Typically 
only half of the 75 eingenvalues of the matrix are finite. I use singular value decomposition methods to perform the 
inversion @]. Once a new p trial is generated, the procedure is iterated until an acceptable x 2 is achieved, which 
I define to be % 2 < 100, more than good enough to make the spectra indistinguishable in practice. Convergence 
typically requires around 20 iterations. In the next section I describe the results of the calculations. 

IV. RESULTS AND CONCLUSIONS 

I choose a target model that is a good fit to the Boomerang and MAXIMA-1 observations, with the following 
parameters: 



total 



l 



1 The structure of the CMBFAST software is particularly well suited to this calculation, as the Boltzmann integration needs 
to be run only once to calculate the entire derivative matrix. 
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n h h 2 = 0.027 

&CDM = 0.3 
fl A = 0.645 
n„ = 
h = 0.7 
n = 1 
t = 0.2 

r = (11) 

Here /i is the Hubble constant in units of lOOkm/s, n is the scalar spectral index, r is the reionization optical depth, 
and r is the tensor/scalar ratio. The procedure is to change one or more of the parameters of the target model to 
produce a "true" model, then reproduce th e Cg spectrum of the target model by modifying the underlying matter 



power spectrum as described in Section HI. This choice of target model has two peculiarities. First, the baryonic 
density is outside the range preferred by Big Bang Nucleosynthesis (BBN), VL h h 2 = 0.019 ± 0.0024 |g||. Second, 
the reionization optical depth is large, r = 0.2. Figure 1 shows the target Ci spectrum and the Ci spectrum resulting 
from setting il^h 2 = 0.019 in accordance with BBN limits. The power spectrum necessary to change the Ci spectrum 
of the true model into one indistinguishable from that of the target model is shown in Figure 2. The perturbation to 
the power spectrum is quite large, and would be difficult to produce via inflation. However, such a power spectrum 
is in principle allowed by existing data. Another choice is to change the reionization optical depth r instead of the 
baryon density. Fig. 3 shows the target model and a model with r = 0. The power spectrum required for the true 
model to mimic the target model is shown in Figure 4. In this case, the power spectrum is dominated by a single 
large feature at long wavelength, a situation that could potentially be realized within an inflationary context. Finally, 
Figures 5 and 6 show the results with f2„ = 0.1, corresponding to a neutrino mass of m„ = 4.5 eV. 

It is clear that allowing for an arbitrary density power spectrum renders attempts at parameter estimation from the 
CMB temperature anisotropy alone a hopeless task. It is necessary to use complementary measurements to constrain 
the density power spectrum. If measurements of the CMB polarization are available in addition to measurements of 
the temperature anisotropy, the task of "fooling" CMB parameter estimation efforts becomes much more difficult. 
Polarization is a tensor quantity, which can be decomposed on the celestial sphere into "electric-type" , or scalar, and 
"magnetic-type" , or pseudoscalar modes. The symmetric, trace-free polarization tensor V a b can be expanded as p2[ 



OO £ 



T ° e=o m -~. 

where the Y (f ' B * , are electric- and magnetic-type tensor spherical harmonics, with parity (— 1)^ and (— 1)^ 



(12) 
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spectively. Unlike a temperature-only map, which is described by the single multipole spectrum of Cj's, a tempera- 
ture/polarization map is described by three spectra 

a Im\ ) = CVy, (\af m \ J=C E e, (\af m \ ) = (13) 
and three correlation functions, 

{ a Jm a fm) = Cce, (aJm a Fm) = C(TB)i, {afm a lm) = C(EB)i- (14) 

Parity requires that the last two correlation functions vanish, C(tb)i = C(EB)i = 0, leaving four spectra: temperature 
Cti-, E-mode Cei, B-mode Cbi , and the cross-correlation Cce- Since scalar density perturbations have no "handed- 
ness," it is impossible for scalar modes to produce B-mode (pseudoscalar) polarization. Only tensor fluctuations (or 
foregrounds p3| ) can produce a B-mode. In the cases I have considered in this paper, calculation of the polarization 
spectra reveals that models with indistinguishable temperature fluctuations have very distinct patterns of polarization, 
so that a sufficiently sensitive measurement of the polarization power spectra (such as that which will be produced by 
the Planck satellite) would reveal any inconsistency in the assumption of a power-law density fluctuation spectrum. 
Figure 7 shows the E-mode polarization spectra for the case considered in Fig. 1, in which the baryon density is 
varied. Fig. 8 shows the cross-correlation spectrum Cti for the same case. The B-mode vanishes. Table I shows the 
X 2 over cosmic variance errors for all three CMB multipole spectra (temperature, E-mode and cross-correlation) for 
the models considered in this paper. In all three cases, although it is possible to mimic the Cti spectrum of the target 
model to within cosmic variance, the polarization spectra produced are poor fits to the respective target polarization 
spectra. Additional information on the density power spectrum can be obtained from large scale structure data. The 
prostpect for independent constraint of the density power spectrum from a combination of SDSS data and CMB 
polarization was considered by Wang et al. 124] . 
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Model 


X 2 : Temperature 


X 2 : E-mode 


X 2 : cross-correlation 


fl b h? = 0.019, Q A = 0.661 


80 


2.2 x 10 b 


4.9 x 10 Y 


T = 


76 


500 


2.1 x 10 4 


Q„ = 0.1, fl A = 0.545 


97 


3.5 x 10 4 


1.8 x 10' 



TABLE I. x 2 f° r various model fits, relative to cosmic variance errors over Ni = 1500 multipole moments. All \ 2 values 
are calculated relative to the target model, with Q, b h 2 — 0.027, Qa = 0.645, and r = 0.2. Qtotal = 1 in all cases. While the 
temperature multipole spectra are excellent fits, the fits of the polarization spectra are poor. 



In this paper I have shown by construction that it is possible to mimic the effect of changes in fundamental 
cosmological parameters on the CMB by changes in the density power spectrum. With a Nk = 75 parameter fit of 
the power spectrum, a CMB multipole spectrum with significantly shifted parameters can be fit to within cosmic 
variance errors, so that even an ideal measurement would be unable to resolve the degeneracy. It is therefore possible 
for estimates of cosmological parameters based on CMB temperature fluctuations alone to be substantially biased 
by features in the density power spectrum. This degeneracy cannot be resolved without access to complementary 
measurements such as observation of CMB polarization or information about large scale structure such as that which 
will be provided by the Sloan Digital Sky Survey. 
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FIG. 1. Ce spectra for target model (black, dashed) and model with Q\>h 2 — 0.019 (red, solid). The spectra differ by 
\ 2 = 1.4 x 10 5 over 1500 degrees of freedom, with cosmic variance errors. The dotted lines show the cosmic variance "envelope" 
for the target model. Error bars are the Boomerang and MAXIMA- 1 measurements, shown for comparison purposes only. 
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FIG. 2. Density power spectrum which mimics the target model {Q.< h 2 = 0.027) with parameters from the true model, 
fib/i 2 = 0.019. With respect to cosmic variance errors, the fit has x 2 = 81 over 1500 degrees of freedom in the Ce spectrum. 
On this plot SP/P = P(k)/k - 1. 
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FIG. 4. Density power spectrum which mimics the C' e s of the target model (r = 0.2) with parameters from the true model, 
r = 0. The final fit has X 2 = 76 over 1500 degrees of freedom. 
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FIG. 6. Density power spectrum which mimics the C' t s of the target model (fi„ = 0) with parameters from the true model, 
Q„ = 0.1. The fit has x 2 = 97 over 1500 degrees of freedom. 
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FIG. 7. E-mode polarization spectra for the target model (Q-^h 2 — 0.027, black, dashed) and the true model (Q,^h 2 = 0.019, 
red, solid), after the power spectrum fit. The dotted lines, barely visible in this plot, show the cosmic variance "envelope" 
of the true model. Even though the temperature multipole spectra of the true and target model are indistinguishable, the 
polarization multipole spectrum differs strongly. 
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FIG. 8. Cross-correlation polarization spectra for the target model (black, dashed) and the true model (Q-^h 2 — 0.019, red, 
solid), after the power spectrum fit, with dotted lines indicating cosmic variance errors. As with the E-mode polarization 
spectrum, the two curves are a poor fit. 
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